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^vq ■ Abstract In the author's article "Algebraic unsolvability of problem of absolute stability 

of desynchronized systems" (Automat. Remote Control 51 (1990), no. 6, pp. 754-759), it 
was shown that in general for linear desynchronized systems there are no algebraic criteria of 
absolute stability. In this paper, a few misprints occurred in the original version of the article 
(f) ■ are corrected, and two figures are added. 

<N . 

^ 1 INTRODUCTION 

. In complex control systems containing sampled-data elements, it is possible that these elements 

^Hh ! operate asynchronously. In some cases asynchronous character of operation of sampled-data ele- 

ments does not influence stability of system. In other cases any small desynchronization of the 
updating moments of sampled-data elements leads to dramatic changes of dynamics of a control 
system, and the system loses stability pp. Last years there is begun (see, e.g., [IHS]) intensive 
studying of the effects connected with asynchronous operation of control systems; both necessary, 
and sufficient stability conditions for various classes of asynchronous systems were obtained. At 
the same time no one succeed in finding general, effectively verified criteria of stability of asyn- 
chronous systems, similar to known for synchronous systems [6]. The problem on stability of linear 
asynchronous systems has appeared more difficult than the problem on stability of synchronous 
systems. In the paper, attempt of formal explanation of complexity of the stability analysis prob- 
lem for linear asynchronous systems is undertaken. It is shown that there are no criteria of absolute 
stability of linear asynchronous systems consisting of a finite number of arithmetic operations. 
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^ I 2 STATEMENT OF PROBLEM 

H ■ 

Consider a discrete-time linear control system whose dynamics is described by the vector differ- 
ence equation 

x(n) = A{n)x{n - 1) (n = 1, 2, . . .), (1) 

where x{n) = {xi(n), X2(n), . . . ,x/v(n)} is the state vector of the system and A(n) = (aij(n)) is a 
square matrix of dimension N with the elements aij(n). 

The system ([1]) will be called synchronous if A(n) = const. If A(n) ^ const, and the set 
{A(n) : n = 1,2,...} consists of finitely many elements Ai,A%, . . . , Am, then the system (pQ) will 
be called asynchronous or desynchronized. 

Let 21 = {Ai, A2, ■ ■ ■ , Am} be a finite totality of square matrices of dimension N. The system 
([T]) will be called absolutely stable with respect to the class of matrices 21 (cf. [7j) if there exists 
c = c(2l) such that for any sequence of matrices A(n) £ 21 the following estimates hold: 

\\A(n)A(n - 1) • • • A(l)z|| < c\\x\\ (n = l,2, ...). (2) 
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Let us call the system (pQ) absolutely exponentially stable with respect to the class of matrices 
21 if there exist c = c(2l) and q = 5 (21) < 1 such that for any sequence of matrices A(n) G 21 the 
following estimates hold: 

\\A(n)A(n- 1) ■ ■ ■ A(l)x\\ < cq n \\x\\ (n = l,2,...). (3) 

If the class of matrices 21 consists of the square matrices A\ = (auj), Ai = [a^ij), ■ ■ ■ , Am = 
{o-Mij) of dimension N then, for its description, it suffices to specify MN 2 numbers: am, 0112, 
. . . , aiNN, 0211, 0212, • • • , a2NN, ■ ■ ■ , clmu, a A/12, • • • , &mnn ■ Therefore, each class 21 consisting 
of M square matrices of dimension N can be treated as a point in some space 9JT(M, N) = R MN . 
Denote by S(M, N) the set of those classes 21 in the space 9Jt(M, N) with respect to which the 
system ([1]) is absolutely stable. By E(M, N) we denote the set of those classes 21 G 9Jt(M, N) with 
respect to which the system ([I]) is absolutely exponentially stable. 

Now, the problem of studying the absolute stability of the system (pQ) can be reformulated as the 
problem of description of the sets S(M,N) and E(M,N); the simpler in some sense the structure 
of the sets S(M,N) or E(M,N) the easier to obtain a criterion of absolute stability or absolute 
exponential stability. 

The sets 5(1, N) and E(l, N) allow a simple description. Indeed, each class 21 G £DT(1, N) consists 
of a single matrix. Therefore, we need to obtain conditions of stability or asymptotic stability of 
some difference equation x(n) = Ax(n — 1). The Routh-Hurwitz stability criterion [6] allows to 
represent these conditions as a finite system of polynomial inequalities including the elements aij 
of the matrix A. Verification of the obtained inequalities can be performed by a finite number of 
arithmetic operations over the elements of the matrix A. In other words, the question whether an 
arbitrary class 21 = {A} belongs to the set S(1,N) or E(1,N) may be resolved by a finite number 
of arithmetic operations. 

Is it possible, for M > 2, by a finite number of arithmetic operations to resolve the question 
whether an arbitrary class 21 G 9Jt(M, N) belongs to the set S(N,N) or E(M,N)7 The answer to 
this question will be given in the next section. 

3 MAIN RESULT 

Let u = {u\, U2, ■ ■ ■ , ul} be an element of the coordinate space R L . A finite sum p(u) = 
Y^Phi2---iL ul i ul 2 ' ' ' w ith numerical coefficients 'Pi x i i ...% Jj is called a polynomial in variable u G R L . 
A set U C R L is said to have the ^-property [8] if there exists a finite number of polynomials 
Pi(u), . . . , Pk{u), pk+i(u), . . . , pk+i{u) such that U coincides with the set of elements u G R L 
satisfying the condition 

pi(u) > 0,... ,pk(u) > 0, Pk+x(u) = ... =pk+i(u) = 0. (4) 

A set U is called semialgebraic [8] if it is a unity of a finite number of the sets possessing the 
5^4-property. 

Theorem 1. Let M , N > 2. If a subset U of the space Wl(M, N) satisfies conditions E(M, N) C 
U C S(M, N) then it is not semialgebraic. 

The proof of the theorem is given in the Appendix. 

Semialgebraicity of a set is equivalent to the existence of a criterion (consisting in verification of 
a finite number of the conditions of the form @) which allows by a finite number of arithmetic op- 
erations of addition, subtraction, multiplication and comparison of numbers to establish belonging 
of an element to a given set. As seen from Theorem [TJ neither the set S(M, N) nor the set E(M, N) 
are semialgebraic. So, the meaning of Theorem [T] is that in general, by a finite number of arithmetic 
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operations, it is impossible to ascertain whether a desynchronized system (pQ) is absolutely stable 
(absolutely exponentially stable) or not. 

The problem on the existence of algebraic criteria of stability is acute also for classes of desyn- 
chronized systems different from those considered above. For example, in the theory of continuous- 
time desynchronized systems there arises the problem of stability of the so-called regular sys- 
tems [HE] (he., the systems with the infinite number of updating moments for each component). 
The discrete-time system ([I]), and the related to it sequence of matrices A(n), will be called regular 
if each matrix Ai from the class 21 = {A\,A 2 , • • • ,Am} appears in the sequence {A(n)} infinitely 
many times. Denote by r(n) the greatest integer r having the property: the set of matrices 
{^4(1), A{2), . . . , A(n)} can be decomposed in r subsets {.A(l), . . . , A(m)}, {A(ni + 1), . . . , A(ri2)}, 
. . . , {A(n r _i + 1), . . . , A(n)} such that each of them contains all the matrices A\,A 2 , ■ ■ ■ ,Am- 
Clearly, the system ([1]) is regular if and only if r(n) — > oo. 

The system (pQ) will be called absolutely exponentially regularly stable with respect to the class 
of matrices 21 = {Ai, A2, . . . , Am} if there exist c = c(2l), q = g(2l) < 1 such that for any regular 
sequence of matrices A(n) E 21 the following estimates hold: 

\\A(n)A(n -!)■■■ A(l)x\\ < cq r{n) \\x\\. 

Denote by R(M, N) the set all the classes 21 E 9K(M, N) with respect to which the system ([I]) is 
absolutely exponentially regularly stable. 

Theorem 2. Let M,N > 2. Then the set R(M,N) is not semialgebraic. 

To prove Theorem[2]it suffices to note that the set R(M, N) contains E(M, N), and is contained 
in S(M, N). Then by Theorem [T] it is not semialgebraic. In other words, for M, N > 2 there are no 
semialgebraic criteria of absolutely exponentially regular stability of discrete-time desynchronized 
systems. 

4 ADDENDUM 

As shown above, the problem of absolute stability of the system (prj can be reduced to the 
analysis of behaviour of infinite products of the matrices A(n) E 21. Theorem [3] below reduces the 
same problem to the descriptive-geometric question on existence in the space R such a norm in 
which each matrix A\, A2, ■ ■ ■ , Am is contractive. 

Theorem 3. The system ([T]) is absolutely stable in a class of matrices 21 = {A±, A2, ■ ■ ■ , Am} 
if and only if there is a norm || • || in R N for which the following inequalities hold: 

||A 1 ||,||A 2 ||,...,p M || < 1. (5) 

The system ([T|) is absolutely exponentially stable in a class of matrices 21 if and only if there is 
a norm || • || in R N and a number q < 1 for which the following inequalities hold: 

\\A 1 \\,\\A 2 \\,...,\\Am\\ < Q- (6) 
The proof of the theorem is given in the Appendix. 

Several important properties of the sets S(M, N) and E(M, N) follow from Theorem [3l For 
example, the set E(M,N) is open in 5DT(M, N); the set E(M,N) belongs to the interior of the set 
S(M, N). 

Due to openness of the set E(M,N), if the system ([T|) is absolutely exponentially stable with 
respect to some class 21° = {A]*,^' • • • >^m} then it is also absolutely exponentially stable with 
respect to any class 21 = {A\,A2, ■ ■ ■ , Am} of matrices Ai sufficiently close to the corresponding 
matrices A® (i = 1, 2, . . . , M). 
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Theorems [T] and [3] imply that the problem of construction, for a given set of square matrices, of 
a norm satisfying conditions © or (|6|) is algebraically unresolvable. 

Theorem [T] states that in general there are no effective criteria of absolute stability of desyn- 
chronized systems (pQ). Nevertheless, such criteria may exist for some particular desynchronized 
systems. Let us present examples. 

Example 1. Denote by 9t(iV) the subset of the space Wl(N, N) consisting of the classes 21 = 
{A±, A2, ■ ■ ■ , An} of matrices A; of the form 
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The problem on absolute stability of the system (JTJ) with respect to the classes 21 £ D\(N) arises 
in [lHS] in the process of study of continuous-time systems with a special types of desynchronization 
of updating moments. 

Theorem 4. The system ([1]) is absolutely stable with respect to the class of matrices 21 = 
{^4i, ^2} G £H(2) of the form ([7]) if and only if one of the following system of relations holds: 

a) an = 1, ai2 = 021 = 0, 022 = 1; 

b) an = 1, ai2 = 0, 022 = — 1, 021 is arbitrary; 

c) an = —1; «12 *s arbitrary, an = 0, 022 = 1/ 

d) an = a 2 2 = — 1, < ai 2 , a 2 i < 4; 

e) I on I < 1, 1 022 1 < 1, -(1 - |an|)(l - I022I) < 012021 < (1 - an)(l - 022)- 

The proof of Theorem U] is cumbersome and so is skipped. Let us point out that the criterion 
of absolute stability of the system ([T]) with respect to the classes of matrices from 9t(2), given in 
Theorem HI is semialgebraic. 

Example 2. Denote by 9\+(N) the subset of the space 9H(iV, N) consisting of the classes 
21 = {A 1} A 2 , ■ ■ ■ ,A N } of matrices A { of the form © for which ay > 0. The set S(N, N) n fR+(N) 
is semialgebraic. The criterion of absolute stability of the system (TjQ) with respect to the classes of 
matrices 21 from 9t+ (N) consists in verification that the maximal eigenvalue of the matrix A = [aij ) 
does not exceed 1. This assertion is proved similarly to Theorem 2 from the first part of pQ. 



APPENDIX 



1. Proof of Theorem [7] suffices to present for the case M = N = 2. The idea of proof is simple. 
We construct two families of matrices depending on the real parameter t £ [-1,1]: 

gw - a - <') ( I ) . *w = (i - ( ) • (A.1) 

The set W of all the classes 21 of the form 21 = = {G(t),H(t)} forms in the space £01(2, 2) 
an algebraic set. Suppose that the set U is semialgebraic. Then the set UC\W is also semialgebraic. 
Therefore, by the theorem of Whitney (see, e.g., [9]) about a finite number of connected components 
of a real algebraic set, a neighbourhood of the class 21(0) in U n W should be either empty or 
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Figure 1: A case when the set U n W has infinitely many components of connectedness 

consisting of a finite number of the connected components. We will show that it has infinitely 
many components of connectedness, see Fig. [TJ So, the set U is not semialgebraic. 

Let us pass to the theorem proof. Denote by | • | the Euclidean norm in R 2 : if x = {£,77} then 
|x| = y^ 2 + rj 2 . Consider two families of matrices: 



Lemma 1. Let ip = vr/(2n + 1). Then (PR n P)((p) = -(cos ip)~ l P{ip). 

Lemma 2. Let p = n/(2n). Then (PR m P)(ip) = \ m , n P((p), where |A m+njn | = \X m , n \, |A m , n | < 

1. 

Both lemmas follow from the equality (PR m P)(f) = cos [( 2 ™+ 1 )y] p(y?) (jn = 0,1,...) whose 
validity is justified by direct calculations, see Fig. [2j 




(A.2) 



PR m Px 




Figure 2: Iterations of a point x under the action of the map PR m P 



Lemma 3. Let ip = ir/(2n), and Bi = P{p) or R>i = R(ip) for 1 < i < m. Then 



B m B m _ x ■■■B x = aR\p)P r (y)R s {ip) 



(A.3) 



where \a\ < 1, integers q and s are non-negative, r = 0, 1. 
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Prove the lemma by induction. For m = 1 the assertion of the lemma is evident; suppose that 
it is valid for m = k — 1 > 1. Then for m = k the matrix A = B m B m ^\ ■ ■ ■ B\ can be represented 
as A = BrnaRvp'R^, where \a\ < 1, r = or f = 1, R = R(<p), P = P(<p). 

If B m = R(ip) then A = aR^ +1 P r R s , and for the matrix A the representation (|A.3P holds in 
which a = a, q = q + 1, r = r, s = s. 

If B m = P((f) and r = then A = aPR g+s , and for the matrix A the representation (|A.3[) holds 
in which a = a, q = 0, r = 1, s = q + s. 

If B m = P(<p) and f = 1 then A = aPR^PR s . Here the factor PR^P according to Lemma [2] 
can be replaced by \q )n P. Then A = a\q )n PR s . Therefore, for the matrix A the representation 
(|A.3|) holds in which a = a\q, n , q = 0, r = 1, s = s. In addition, \a\ < \a\ ■ \Xq, n \ < 1 since \a\ < 1, 

\^q,n\ < 1- 

The inductive step is completed. Lemma [3] is proved. 
Corollary. \B m B m ^ ■ ■ ■ B x \ < \P (£) \ . 

The proof of the corollary immediates from the representation (|A.3P and unitarity of the rotation 
matrix R(^p). 

Lemma 4. Let t n = sin ^. Then 2l(t„) G E(2, 2). 

Proof. Let {A(k)} be a sequence of matrices from 2l(t n ). Then, for each k, one of two equalities 
A(k) = G(t n ) or A(k) = H(t n ) holds. By G(t n ) = #(t n ) = /i^(^), where 

/i n = 1 — (sin . Therefore, the product of matrices ^4(1), ^4(2), . . . , A(k) can be represented in 
the form: A(k)A{k - 1) • • • A(l) = ^ k n B k B k ^ ■ ■ ■ B u where B< = P (£) or B t = R (^). Then, by 
Corollary from LemmaH \A(k)A(k — 1) • • • ^4(1)| < /ijj \P (^) | which implies absolute stability of 
the class of matrices 2l(i n )- Lemma H] is proved. 

Lemma 5. Let s n = sin 2 n+i • Then 2l(s n ) S(2, 2) for all sufficiently large n. 

Proof. Clearly the lemma will be proved if, for each sufficiently large n, there can be found a 
sequence of matrices A(k) G d{s n ) such that 



\A(ki)A(ki-l)...A(l)\ -+oo 



(A.4) 



for some fci — > oo. 

Define the sequence of matrices A(A;) as follows: j4[(n + 2)z] = G(s n ), A[(n + 2)i + 1] = . . 
A[(n + 2)i + n] = H(s n ), A[(n + 2)i + n + 1] = G(s n ). Let us set fcj = (n + 2)i + n + 1. Then 



- 1) • • • A(l) = [G(s n )H n (s n )G(s n )} 
Since G(s n ) = v n P ( ^pi) and H(s n ) = u n Q (^ppr), where v n = 1 - 
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2n+l 



, then 



-!)■■• A(l) 



K +2 p I " 

1 2ra + 1 



g n 



7T 



2n + 1 



P 



7T 



2n + 1 



Consequently, by Lemma Q] A(fcj)A(A;j — !)••• A(l) 



P f 2n+i J * s a P r °j ec tor and so 



\^2n+l 



< +2 / cos ) p ( 2^t ) • Reca11 that 



> 1. Therefore, 



\A(ki)A(ki-l) — A(l)\ > 



,n+2 



COS 



2n+l 



(A.5) 
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Direct calculations show that ^ +2 / cos 2 n+i = 1+ 2(2n+i)^ +°( n 2 )- Henceforth, for sufficiently 
large values of n the inequality f ™ +2 / cos > 1 holds. From this and from (1A.5|) the relation 
()A.4p follows. Lemma [5] is proved. 

Let us complete the proof of the theorem. Since by Lemmas H] and [5] 2l(t n ) G £'(2,2), 2l(s n ) G" 
5(2, 2) and E{M, N) C [7 C 5(M, iV) then 2l(i n ) G 17, 2l(s n ) [7. But because of the points t n and 
s n interleave between each other then the set U(~)W contains infinitely many different components of 
connectedness (different classes 2t(i n ) belong to different components of connectedness). Therefore, 
the set U fl W by the theorem of Whitney [U] is not semialgebraic. But since W is an algebraic set 
then the set U is not semialgebraic. Theorem [1] is proved. 

2. Proof of Theorem^ In one side the assertion of Theorem [3] is obvious: absolute stability and 
absolute exponential stability of the system (pQ) with respect to the class 21 = {A\,A 2 , ■ ■ ■ , Am} 
immediately follow from inequalities © and ([6]). 

Let us show that absolute exponential stability of the system (pQ) with respect to the class of 
matrices 21 = {A\,A 2 , . . . ,Am} implies ©. Let, for some q < 1, the relation §3§ be valid. Set 
kq(x) = \x\, >c n {x) = q~ n max \B\B 2 ■ ■ ■ B n x\ (n > 1), where | • | is the Euclidean norm on R N , and 
the maximum is taken over all possible collections of the matrices B\, B 2 , . . . , B n G 21. Define the 
norm || • || as follows: ||a;|| = sup n>0 >c n (x). 

The function ||x|| is semiadditive and due to ([3]) it satisfies the relations |x| < ||x|| < max{l, c} \ x\. 
Henceforth, ||x|| = if and only if x = 0. Other properties of a norm are obvious for || • ||. 

Let us justify inequalities ([6]). Clearly, for any k = 1,2, ... , M and n = 0, 1, . . . the estimate 
x n (>lfcx) < qx n+ i(x) is valid. Therefore, 

\\A k x\\ = supx n (A k x) < gsupx n+ i(x) = supx n (a;) < q\\x\\. 

n>0 n>0 n>l 

From here 1 1 ^4.^ 1 1 < q for k = 1,2, , , . ,M. Inequalities ([6]) are proved. 

Construction of the norm || • || and the proof of inequalities © in the case of absolute stability 
of the system ([T]) are carried out similarly. Theorem [3] is proved. 
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POST SCRIPTUM 

In the foregoing text, a few misprints in the proof of Theorem Q] occurred in the original journal 
version of the article |10| were corrected, and two figures were added. The improved text was 
included in the monograph |11| . Generalization of the presented results can be found in [12J. 
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